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I.  INTRODUCTION 


During  May  -  June  1083,  the  Ballistic  Research  Laboratory  (BRL)  and  the  Human 
Engineering  Laboratory  (HEL)  conducted  an  experiment  at  the  joint  BRL-HEL 
Command  Post  Exercise  Research  Facility  (CPXRF).  The  overall  purpose  of  the 
experiment  was  to  demonstrate  the  feasibility  of  utilizing  the  automated  techniques  of 
the  facility  to  run  fire  support  control  experiments.  For  the  first  time,  military  players 
interfaced  with  the  newly  developed  Artillery  Control  Environment  (ACE)  software. 

To  demonstrate  this  capability,  the  effects  of  message  intensity  and  communi''ation 
degradation  on  the  FIST  HQ  ability  to  perform  fire  support  coordination  were  studied. 
One  of  the  measures  of  performance  used  to  quantitatively  evaluate  these  effects  was  the 
distribution  of  message  service  times  for  the  fire  support  team.  Service  time  has  been 
defined  to  be  the  time  required  by  the  FIST  HQ  to  service  a  message  beginning  at  the 
time  an  acknowledgement  (ACK)  is  sent  from  the  fire  support  team  digital  message 
device  (FIST  DMD),  indicating  receipt  of  a  message,  until  the  time  a  response  message  is 
first  transmitted.  This  measure  indicates  the  time  a  message  spends  in  the  FIST  DMD 
message  queue  combined  with  the  processing  and  decision  time  of  the  FIST  HQ. 

An  analysis  of  variance  (ANOVA)  procedure  performed  on  the  FIST  HQ  message 
service  time  data  indicated  that  level  of  intensity,  communication  degradation  level, 
message  type,  team,  and  replicate  (of  treatment  combination)  all  had  a  significant  effect 
on  the  FIST  HQ  service  time.*  Only  the  service  time  distributions  of  fire  request  (FR) 
and  artillery  target  intelligence  (ATI)  messages  were  modeled  and.  because  message  type 
was  significant,  both  types  were  modeled  separately.  These  were  the  only  message 
types  (other  than  freetext  messages  which  were  excluded)  that  were  not  automatically 
forwarded  by  the  FIST  DMD  while  it  was  operating  in  the  automatic  mission  mode.' 
The  service  times  for  fire  request  messages  were  modeled  based  on  the  data  collected  for 
the  36  two-hour  cells  (all  significant  factors  considered).  While  there  was  marked 
variation  among  the  median  ATI  service  times  by  experimental  condition,  ihe  service 
times  for  ATIs  were  combined  over  team  and  replicate  and  modelr  !  based  on  9  eight- 
hour  cells.  This  was  expedient  due  to  a  short  time  frame. 

The  model,  or  models,  describing  the  distribution  of  FIST  HQ  message  service 
times  can  eventually  be  incorporated  into  large  simulation  codes  of  Field  Artillery- 
command,  control,  and  communications  (C^)  and  a  FIST  HQ  simulator  which  is 
currently  under  development. 


^  Smith,  Jill,  GiTBOvicki,  Jock  O  ,  tt  al ,  *Fire  Sappon  Team  Expenmeit,'  BRI/-MR-3423,  December  1984.  AD  AISO  297 

After  the  FIST  HQ  had  received  aa  laitial  &re  reqoeet  and  decided  who  would  handle  the  mij.Mon.  the  Gre  mission  was 
forwarded  ia  the  antomatie  mitaioa  mode  All  aobteqoeat  messa(ei  for  that  mission  were  antomatically  routed  through  the  FIST 
DMD  If  a  message  was  not  acknowledged  in  fonr  transmissions,  the  FIST  DMD  operator  was  notiGed,  the  message  was  placed  in 
his  menage  qnene  and  then  transmitted  mannally 


U.  MODELING  APPROACH 

IdentificatioD  of  an  appropriate  functional  form  with  which  to  model  the  FIST  HQ 
message  service  time  data  involved  at  least  three  major  steps.  First,  was  the  selection 
of  a  preliminary  model  which  was  based  on:  1)  an  understanding  of  the  origin  of  the 
data  to  help  choose  the  family  of  functional  forms  to  be  tested,  and  2)  graphical 
representation  of  the  data  and  observation  of  specific  functions  that  would  produce 
forms  similar  to  the  ones  exhibited  by  the  graphs.  After  the  selection  of  a  preliminary 
model,  the  second  step  was  to  determine  the  relevant  model  parameter  values  which,  in 
this  case,  had  to  be  estimated  from  the  empirical  data.  The  third,  and  final,  step  in 
the  model  identification  process  was  the  verification  of  the  model.  This  was 
accomplished  by  using  appropriate  goodness-of-fit  tests  to  compare  the  observed 
distributions  of  ATI  and  fire  request  messages  with  the  hypothesized,  or  theoretical, 
model.  For  ATI  messages,  the  chi-square  (x^)  goodness-of-fit  test  was  used  for 
comparison  of  the  observed  frequency  distribution  under  varying  conditions  with  the 
proposed  theoretical  frequency  distribution.  For  fire  request  messages,  due  to  the  much 
smaller  sample  sizes  and  the  non-normality  exhibited  by  the  observed  data,  a  specialized 
Kolmogorov-Smirnov  test  was  usee  ‘  o  test  the  fit  of  the  same  model. 

A.  Selection  of  the  Preliminary  Model 

One  of  the  fundamental  problems  of  modeling  any  type  of  data  is  the 
determination  of  the  criteria  to  select  an  appropriate  family  of  densities  from  which  to 
identify  a  reasonable  functional  form  to  fit  to  the  observed  data.  Typically,  service  times 
exhibit  considerable  variation  and  may  usually  be  regarded  as  random  variables.  It  is 
therefore  natural  to  represent  service  times  in  terms  of  a  random  process.  The 
exponential  family  of  densities,  frequently  invoked  in  modeling  studies  to  describe 
service  time  processes,  was  selected  as  the  preliminary  functional  form  for  this  model. 
This  was  primarily  based  upon  the  descriptive  properties  exhibited  by  the  FIST  HQ 
observed  service  time  frequency  distributions. 

The  observed  frequency  distributions  of  the  two  message  types,  fire  request  and 
artillery  target  intelligence,  have  been  graphically  represented  as  the  histograms  in  the 
figures  of  Appendices  A  and  B,  respectively.  The  histograms  for  either  message  type 
do  not  exhibit  the  property  of  symmetry  but,  in  general,  taper  off  to  the  right  of  the 
peak  of  the  distribution  and  can  be  described  as  asymmetric. 

There  is  a  very  distinct  "rise  and  fall”  pattern  evident  in  all  the  artillery  target 
intelligence  message  histograms;  this  same  pattern  is  repeated  throughout  most  of  the 
histograms  of  the  fire  request  messages.  However,  the  shape  of  this  "rise  and  fall" 
pattern  of  the  observed  frequency  distributions  varied  not  only  between,  but  also  within, 
the  two  message  types.  This  shape  is  a  descriptive  measure  of  the  amount  of  dispersion 
(variation)  within  the  observed  data  and  is  indicative  of  how  strongly  each  distribution 

It  shonid  be  noted  for  the  purpose  of  comparison  that  the  x-nxis  for  FRs  and  ATIs  is  scaled  the  same,  0.0  seconds  to  50  0 
seconds,  while  the  FR  y-axis  (0  to  6)  is  scaled  at  one-half  the  ATI  y-axis  (0  to  12)  Also,  for  both  FRs  and  ATIs,  whenever 
observations  occurred  at  >  50  0  seconds,  they  were  grouped  into  a  single  cell  and  graphed  at  x=50  0  seconds.  The  actual  range  of 
the  data  within  this  cell  has  been  annotated  at  this  point 
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concentrated  about  some  particular  central  value,  in  this  case,  the  median  service  time. 
The  measure  of  dispersion  should  be  large  when  the  spread  of  the  distribution  about  the 
median  is  large  and  should  be  small  when  the  spread  is  negligible.  Variations  in  shape 
among  the  distributions  were  not  completely  unexpected  since  there  was  marked 
heterogeneity  of  variance  among  the  two-hour  cells  of  the  original  observations.' 

Although  not  evident  from  the  histograms,  the  origin  of  each  observed  frequency 
distribution  was  translated  on  the  X-axis  from  an  unknown  location  point  (threshold)  to 
zero.  This  location  point  is  defined  to  be  the  smallest  observed  service  time  for  a  set  of 
responses  under  particular  experimental  conditions.  The  value  of  the  location  point 
ranged  between  3.0  and  6.0  seconds  for  ATIs  with  the  median  and  mode  being  4.0 
seconds,  while  for  fire  requests  (FR)  it  ranged  between  4.0  and  25.0  seconds  with  8.0 
seconds  as  both  the  median  and  mode. 

The  probability  density  function  selected  to  initially  model  the  FIST  HQ  service 
times  had  to  be  able  to  profile  the  descriptive  properties  mentioned.  In  summary,  the 
distribution  had  to  satisfy  the  following  criteria: 

1)  be  asymmetrical, 

2)  include  a  parameter,  or  parameters,  to  measure  both  the  dispersion  and 
the  rate  at  which  the  "tair  of  the  function  tapered  off, 

3)  include  a  measure  of  location. 

The  three-parameter  WeibuII  distribution  from  the  exponential  family  of  densities 
satisfied  these  criteria.  If  the  random  variable  service  time,  call  it  X,  has  a  three- 
parameter  Weibull  distribution,  then  its  probability  density  function  fx  (x)  of  X  has  the 
form 


if  X  >  1/ 
if  X  <  1/ . 

The  three  constants  /?  >  0,  a  >  0,  and  r/  >  0  are  the  parameters  of  the  distribution. 
The  parameter  /?  determines  the  shape  of  the  density,  1/a  is  a  scale  parameter  that 
measures  the  steepness  of  the  fall  of  the  distribution,  and  the  parameter  v  provides 
information  about  the  smallest  possible  value  of  the  random  variable  X  (service  time). 

The  roles  of  /3  and  o  are  illustrated  in  Figures  1  and  2.  In  Figure  1  fx(x)  is  graphed 
for  u  =  0.0,  /?  =  1.5,  Q  =  0.5,  2.0,  and  in  Figure  2  fx(x)  is  graphed  for  v  =  0.0,  = 

5.0,  a  =  0.5,  2.0. 


The  Weibull  Density  Function  for 


Parameter  Estimation  for  the  WeibiiH  Distribution 

Once  the  preliminary  model  for  the  FIST  HQ  service  time  data  had  been  selected, 
the  next  step  was  to  determine  the  relevant  parameter  values.  Recall  that  the 
parameters  I3,  a,  and  t/  of  the  Weibull  distribution  denote  shape,  scale  and  location, 
respectively. 

Analytical  estimates  of  the  three  distribution  parameters,  when  all  are  unknown, 
should  be  obtained  iteratively  for  which  the  method  of  maximum  likelihood  is 
appropriate.  Since  the  maximum  likelihood  estimators  cannot  be  obtained  explicitly, 
another  method  proposed  by  Dubey^’*  provides  a  simple  point  estimator  for  i/  and  also 
allows  one  to  explicitly  obtain  percentile  estimators  for  I3  and  a.  These  estimates  may 
be  considered  satisfactory  approximations  of  a,  and  u  to  be  used  as  initial  guesses  in 
the  solution  of  the  likelihood  equations. 

Suppose  that  X(,j,  x^2],  '  ‘  ‘  i  ^(n)  represent  n  ordered  observations  of  a  random 
sample  of  the  FIST  HQ  service  times  from  the  Weibull  population  with  fx(x)  given  by 
Eq.  (1).  Then  Dubey  proposed  that  a  simple  estimator  for  u  is  of  the  form 

X(t)  X(fc)  -  X^(j) 

X(1)  +  X(^)  -  2X(j)  ’ 

where  Xjj)  and  Xj^j  are  any  of  the  ordered  observations  such  that  xji)  <  xjjj  <  (Xjij  Xj^jp. 
Throughout  the  process  of  computing  the  point  estimates  of  u  for  each  cell  modeled,  Xjjj 
and  X(|j)  are  selected  as  the  ordered  observations  corresponding  to  the  25th  and  75th 
sample  percentiles,  respectively. 


Once  1/  has  been  estimated,  a  percentile  ^timator  of  the  shape  parameter,  P,  can  be 
explicitly  obtained  based  on  any  two  sample  percentiles.  However,  Dubey®  showed  that 
the  17th  and  97th  sample  percentiles  asymptotically  yield  the  best  percentile  estimator 
of  0.  These  are  the  sample  percentiles  used  in  the  estimation  of  the  shape  parameter. 
The  location  parameter  v  is  neither  known  to  be  0  nor  has  it  yet  been  set  equal  to  0  by 
an  appropriate  transformation  of  the  data.  A  slight  modification  of  the  percentile 
estimator  of  13,  based  on  two  ordered  sample  observations  from  the  Weibull  population 
and  assuming  t/  is  known  but  different  from  0,  is  given  by; 

.  ^  In  (-  ln(l-pi))~ln  (-  In  (1 -pg)) 

In  (yp,  -  i')  -  In  (yp^  -  u) 

where 

0  •<  P|  ■<  P2  1;  Pi  ~  .1673,  P2  =  .9737, 
yp,  =  yi,  i  =  1,2 

and 

“2 - 

Dubfy,  Satya  D  ,  "Hypsr-ffficifnt  Eatimator  of  the  Loratioa  Parameter  of  the  Weiboll  Laws,”  Naval  Researrh  l.oeistirs 
Quarterly.  Vol  13,  1966 
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Dubey,  Satya  D  ,  ’Some  Percentile  Estimators  for  Weibull  Paremeten,”  Terhuoitietiics  Vol.  9,  1967 


Xp,  the  100  p  per  cent  percentile  of  the  sample,  is  defined  as 

{Xjnpj  if  np  is  an  integer 
X(jnpi+i)  if  is  an  integer 

where  [np]  denotes  the  largest  integer  not  exceeding  np  and,  as  previously  stated,  n  is 
the  size  of  the  sample. 

Finally,  when  it  is  assumed  that  the  shape  parameter  is  known  and,  again,  that  i/  is 
known  but  different  from  0,  then  a  modified  single  observation  best  percentile  estimator 
of  a  is 

a  =  y  /  In*/^  (1  -  p)-‘ ,  (6) 


where  y  is  any  sample  percentile  corresponding  to  the  cumulative  probability  p  (0  <  p 
<  1).  For  the  sake  of  consistency,  p  was  set  equal  to  pj  and  y  set  equal  to  yp^,  where  P| 
and  ypi  were  previously  defined  for  the  development  of 

The  point  estimates  of  d,  and  i/,  calculated  using  the  above  procedures,  are  given 
in  Tables  C-lA  thru  C-3A,  and  C-lOA  for  fire  requests  and  Alls,  respectively.  These 
estimates  were  input  as  initial  guesses  to  the  Biomedical  Computer  Programs  P-Series 
(BMDP)  maximum  likelihood  estimation  program  with  one  exception.  The  point 
estimates  of  the  location  parameter,  if,  for  both  FRs  and  ATIs  closely  approximated  the 
minimum  service  time  observed  in  each  of  the  cells  being  modeled.  Since  the  location 
parameter  has  been  defined  as  the  smallest  order  statistic  of  an  ordered  set  of 
observations,  and  no  prior  knowledge  existed  to  support  the  theory  that  u  should 
assume  some  particular  value  !/„,  it  was  decided  to  fix  each  estimate  of  as  the  smallest 
observed  service  time  within  each  cell.  These  fixed  estimates  of  the  location  parameter, 
along  with  the  maximum  likelihood  estimates  of  ^  and  d,  are  given  in  Tables  C-lB  thru 
C-3B  for  fire  requests  and  in  Table  C-lOB  for  artillery  target  intelligence  messages. 

1.  Fire  Request  Messages 

The  maximum  likelihood  estimates  of  both  0  and  d  for  fire  request  messages  (see 
Tables  C-lB  thru  C-3B)  proved  to  be  quite  varied  within  each  of  the  significant 
experimental  conditions:  level  of  communication  degradation,  replicate,  team,  and  level 
of  intensity.  Large  values  of  8.844  and  35.391  for  0  and  d,  respectively,  occurred  at 
30^  communication  degradation,  replicate  1,  team  3,  and  low  intensity.  A  possible 
explanation  of  why  such  large  estimates  occurred  during  low  intensity  will  be  given 
during  the  description  of  the  parameter  estimates  of  if. 

With  the  value  of  35.391  excluded,  the  maximum  likelihood  estimates  of  d 
appeared  to  be  fairly  consistent  across  the  factors,  ranging  from  a  low  of  approximately 
1.000  to  a  high  of  less  than  15.000.  It  should  be  noticed  in  the  tables  that  small  values 
of  d  and  li  are  usually  paired  and  if  d  is  large,  then  is  also  large.  Recall  that  these  two 
parameters  affect  the  shape  of  the  distribution  (i.e.,  how  highly  concentrated  the 
distribution  is  about  a  selected  measure  of  central  tendency).  The  variance  of  the 


random  variable,  FIST  HQ  fire  request  message  service  time,  provides  information  about 
the  "scatter”  or  amount  of  concentration  of  the  observations.  Since  the  variance  of  a 
WeibuU  distribution  is  dependent  upon  the  two  parameters  a  and  the  wide  range  of 
values  for  the  estimates  of  a  and  ^  is  consistent  with  the  non*homogeneity  of  variance 
of  the  observed  service  times  mentioned  earlier. 

The  parameter  u  serves  only  as  a  measure  of  location  and  does  not  infiuence  the 
shape  of  the  distribution  in  any  way.  Tables  C-lB  thru  C-3B  give  the  maximum 
likelihood  estimates  of  u  for  the  fire  request  service  time  data.  Of  the  four  factors 
involved,  the  biggest  difference  among  the  estimates  of  u  occurred  between  the  two 
teams.  Estimates  of  u  ranged  from  8.000  to  25.000  seconds  for  team  3  and  4.000  and 
12.000  seconds  for  team  4.  There  were  two  extreme  values  for  if  of  25.000  seconds 
during  Rep  (replicate)  I  and  16.000  seconds  during  Rep  II  for  team  3  while  operating 
under  low  intensity  and  30%  communication  degradation.  It  should  be  pointed  out  that 
as  the  level  of  intensity  increased  from  low  to  high,  if  slightly  decreased.  At  low 
intensity  if  ranged  between  5.000  and  25.000  seconds,  at  medium  intensity  it  ranged 
from  5.000  to  14.000  seconds,  and  at  high  intensity  if  ranged  from  4.000  to  12.000 
seconds.  This  may  appear  to  be  an  anomolous  situation,  but  often  a  server  (in  this  case 
the  FIST  HQ)  increases  his  effort  when  under  pressure,  perhaps  due  to  a  long  message 
queue. 

With  the  two  large  estimates  of  if  for  team  3  removed,  the  mean  of  if  with  all 
factors  considered  is  8.412  seconds.  The  mean  if  over  all  associated  leveb  of  each 
experimental  condition  proved  to  be  quite  stable  ranging  from  8.361  seconds  (over  levels 
of  intensity)  to  a  high  of  only  8.524  seconds  (across  teams).  These  estimates  may 
provide  some  insight,  particularly  for  modeling  purposes,  as  to  the  minimum  amount  of 
time  required  by  the  FIST  HQ  to  review  and  transmit  a  fire  request  message;  thb  time 
would  also  include  the  time  the  message  remained  in  the  FIST  DMD  message  queue. 
They  may  help  in  laying  a  theoretical  foundation  upon  which  i/  can  be  fixed  at  some 
particular  iff,. 


Maximum  likelihood  estimates  of  ^  for  ATIs  (see  Table  C-lOB)  ranged  between 
1.001  and  2.075  with  8  of  the  0  estimates  either  greater  than  or  equal  to  1.000  but  less 
than  2.000.  When  is  estimated  to  be  very  close  to  1.000,  the  Weibull  distribution  has 
no  mode  and  the  probability  density  function  decreases  monotically  as  the  random 
variable  service  time  increases.  Where  ^  is  distinctly  greater  than  1.000,  such  as  for  low 
intensity  and  00%  communication  degradation,  the  distribution  is  unimodal. 

With  the  exception  of  the  extremely  large  value  of  d,  12.555,  at  high  intensity  and 
a  30%  communication  degradation  level,  the  d  estimates  exhibited  only  a  slight  degree 
of  variation.  Excluding  the  d  value  of  12.555,  d  ranged  between  4.837  and  0.452.  The 
large  estimate  of  d  at  high  intensity  and  30%  communication  degradation  is  not  really 
surprising  as  it  is  indicative  of  a  much  longer  tailed  distribution  (larger  service  times) 
which  is  not  as  concentrated  about  the  large  median  service  time  of  10.5  seconds  for  this 
cell.*  This  is  consistent  with  the  fact  that  when  the  FIST  HQ  is  operating  under  such 


severe  conditions,  FRs  have  a  greater  priority  and  will  be  processed  before  any  incoming 
ATb  which  are  then  forced  to  remain  in  the  message  queue,  resulting  in  a  more 
widespread  distribution  of  service  times  for  this  particular  cell  than  any  other.  This  is 
evident  from  its  histogram  of  observed  service  times  (see  Appendix  B)  which  range 
between  0  and  93  seconds,  a  range  which  is  almost  half  again  as  large  as  the  next  largest 
range  of  service  times  occurring  at  low  intensity  and  30%  communication  degradation. 

Recall  that  u  was  fixed  as  the  shortest  observed  time  in  which  a  message  could  be 
processed  and  transmitted.  For  ATIs  there  was  very  little  variation  among  the 
estimates  either  within  the  levels  of  intensity  or  communication  degradation;  overall,  i/ 
varied  between  a  low  of  3.000  seconds  and  a  high  of  only  6.000  seconds.  The  largest 
observed  value  of  u,  6.000  seconds,  occurred  at  high  intensity  and  30%  communication 
degradation,  probably  for  the  same  reason  as  mentioned  in  the  previous  paragraph. 

C.  Transformation  of  the  Weibull  Distribution  to  the  Exponential  Distribution 

The  estimates  '0,  a,  and  P,  of  the  three-parameter  Weibull  distribution  provide  a 
reasonable  analytical  interpretation  of  the  effects  that  various  combinations  of  level  of 
intensity  and  level  of  communication  degradation  had  on  the  distribution  of  the  FIST 
HQ  service  times  for  fire  request  and  artillery  target  intelligence  messages.  However, 
considerable  effort  is  required  to  compute  these  estimates  regardless  of  the  procedure 
utilized. 

The  exponential  distribution  was  initially  ruled  out  because  this  distribution  is 
limited  by  the  assumption  of  a  constant  service  rate  (the  function  is  monotonically 
decreasing)  while  the  Weibull  distribution  can  be  written  to  include  variable  (increasing 
and  decreasing)  service  rates.  However,  an  exponential  model  can  be  used  in  situations 
where  the  service  rate  is  not  constant  as  long  as  an  appropriate  transformation  of  the 
Weibull  distributed  random  variate,  FIST  HQ  service  time,  is  made  (see  Appendix  D). 

To  make  such  a  transformation,  the  following  assumptions  about  the  three 
parameters  of  the  Weibull  distribution  bad  to  be  made;  1)  the  parameters  and  r/  are 
known  while  a  is  assumed  unknown,  2)  0  =  1.000,  and  3)  i/  ==  0.000.  In  Section  II-B,  it 
was  shown  that  0  was  equal  to  or  "near”  1.000  about  90%  of  the  time  for  ATIs,  and 
estimated  at  or  "near”  1.000  about  67%  of  the  time  for  FRs.  It  therefore  seemed 
reasonable  to  conclude  that  the  shape  parameter  0  was  indeed  known  and  could  be  set 
equal  to  1.000.  Based  on  this  conclusion  an  exponential  probability  density  function 
with  both  scale  and  location  parameters  was  derived.  To  obtain  the  final,  single 
parameter  exponential  model  each  observed  service  time  within  a  particular  cell  was 
adjusted  by  an  amount  equivalent  to  Xn,i„  =  u  for  that  particular  cell.  Values  of 
X,„5„  =  u  for  the  cells  are  given  in  Tables  C-lB,  C-2B,  and  C'~3B  for  FRs  and  Table  C- 
lOB  for  ATIs.*  This  translation  effectively  shifted  the  origin  of  the  observed  service  time 
distributions  to  X_:„  =  u  =  0.000. 


t'  -  ■  ■ 

Note  tkit  theae  are  the  lame  Taloea  of  that  were  fixed  denog  the  computation  of  maximum  likelihood  r'timatrs  of  O  and  fJ 
for  the  Weibull  diatribntion. 


Since  the  assumptions  have  been  met,  the  three  parameter  Weibull  distribution  is 
equivalent  to  an  exponential  distribution  with  scale  parameter  B  =  1/a.  The 
probability  density  function  of  the  new  random  variate  Y,  adjusted  FIST  HQ  service 
time,  becomes 

B  e-'y,  if  y  >  0, 

0,  if  y  <  0. 

Maximum  likelihood  estimates  of  B  for  varying  communications  degradation  are 
presented  in  Tables  C-4  thru  C-6  and  C-11  (Appendix  C)  for  fire  request  messages  and 
ATIs,  respectively.  Since  the  scale  parameter  a  of  the  Weibull  distribution  was  assumed 

unknown,  the  maximum  likelihood  estimates  of  B  were  not  computed  as  l/a  but  rather 
N  N  _ 

as  ( 1.0/  (  Y]  yi/N)  1  ,  where  y;  =  x;  -  i/  and  Xi/N)  is  the  mean  of  a  cell  of  adjusted 
i=.i  i-l 

service  time  observations.  There  is,  however,  reasonable  agreement  between  the 
maximum  likelihood  estimates  of  0  as  generated  by  the  above  method  and  likelihood 
estimates  of  B  if  they  had,  in  fact,  been  computed  as  1/a.  As  an  example,  in  Table  C- 
lOB  for  ATIs,  with  00%  communication  degradation  and  low  intensity,  the  maximum 
likelihood  estimate  of  a  is  6.582  making  the  estimate  of  ^  =  1/a  =  .152.  Table  C-11 
gives  the  maximum  likelihood  estimate  of  B,  under  the  same  conditions,  as  .150. 

Exponential  fits  to  the  observed  data  sets  using  the  maximum  likelihood  estimates 
of  B  with  V  fixed  at  0.000  are  presented  in  Appendices  A  and  B  as  the  heavily  shaded 
curves. 

D.  Model  Verification 

The  final  step  in  the  model  identification  process  was  the  veriiication  of  the  single 
parameter  exponential  model.  This  was  accomplished  by  performing  an  appropriate 
goodness-of-fit  test  on  the  observed  data  for  each  message  type  to  compare  its 
hypothesized  probability  distribution  with  the  frequency  distribution.  The  null 
hypothesis  for  the  goodness-of-fit  test  is  a  statement  about  the  identity  of  the 
probability  distribution  that  fits  the  FIST  HQ  service  time  data.  Normally,  this 
hypothesized  distribution  is  completely  specified,  including  all  parameter  values. 
However,  as  previously  mentioned,  the  unknown  parameter  of  the  exponential 
distribution  had  to  be  estimated  from  the  data  before  any  goodness-of-fit  test  could  be 
carried  out.  It  should  be  understood  that  rejection  of  the  null  hypothesis  for  any 
particular  case  gives  no  information  about  what  the  true  service  time  population  form 
is,  only  what  it  is  not.  Thus,  the  goodness-of-fit  tests  were  utilized  to  obtain  statistical 
support  for  the  theory  that  the  single  parameter  exponential  distribution  fits  both  the 
FR  and  ATI  message  service  time  data  by  acceptance  of  the  null  hypothesis. 

Two  well-known  statistical  tests  for  goodness-of-fit  were  used  for  the  model 
verification  process.  An  adaptation  of  the  Kolmogorov-Smirnov  (K-S)  test  statistic 
was  used  to  verify  the  model  for  Are  request  messages,  while  the  chi-square  test  was 
used  to  verify  the  same  model  for  artillery  target  intelligence  messages.  Both  tests 
are  based  on  a  comparison  between  the  distribution  of  the  observed  sample  data  and  the 
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theoretical  distribution  that  has  been  hypothesized;  however,  the  two  tests  use  a 
different  basis  for  comparison. 

The  technique  used  for  testing  the  good ness-of- fit  between  the  observed  sets  of 
FIST  HQ  service  times  for  fire  request  messages  and  the  hypothesized  exponential  model 
was  the  Kolmogorov-Smirnov  test.  Its  measure  of  incompatibility  is  based  on  the 
vertical  deviations  between  the  observed  and  the  hypothesized  cumulative  distribution 
functions.  The  null  hypothesis  of  the  K-S  test  specifies  that  for  the  observed  random 
variable,  FR  service  time,  the  cumulative  probability  function  is  some  Fo(x).  This  is 
tested  against  the  alternative  hypothesis  that  the  cumulative  probability  function  of  the 
observed  random  variable  is  different  from  Fo(x).  This  null  hypothesis  was  tested  at  a 
level  of  significance  of  0.05.  One  important  assumption  of  the  K-S  test  is  that  the  sets 
of  observations  come  from  a  completely  specified  continuous  distribution.  However,  the 
scale  parameter  0  of  the  exponential  distribution  was  not  known  and  had  to  be 
estimated  from  the  sample  data,  hence,  the  standard  K-S  test  no  longer  applies  because 
the  commonly  tabulated  critical  values  of  the  test  statistic  become  conservative 
(although  exactly  how  conservative  is  generally  unknown).  In  this  case,  tables  listing 
critical  values  for  a  specialized  Kolmogorov-Smirnov  test,  for  the  exponential 
distribution  when  its  population  mean  is  unknown,  had  to  be  used  in  place  of  the 
standard  tables.^  This  specialized  test  provides  a  good  ness-of- fit  test  which  can  be  used 
with  sample  sizes  that  are  considered  too  small  for  use  with  the  chi-square  test. 

The  Kolmogorov-Smirnov  goodness-of-fit  test  statistic  measures  agreement  as  the 
absolute  value  of  the  largest  vertical  difference  between  the  graphs  of  the  two 
cumulative  relative  frequency  distributions.  It  is  defined  as  D  =  max^lFo(x)  -S(x)l, 
where  S(x)  is  the  sample  (observed)  cumulative  distribution  function,  Fo(x)  is  the 
hypothesized  cumulative  exponential  distribution  function  with  parameter  0,  and  the 
max  over  x  is  the  largest  of  the  differences  |  F0(x)  -  S(x)  |  in  the  neighborhood  of  each 
observed  value  of  x.  Since  the  calculated  value  of  Z)  is  a  measure  of  agreement  between 
the  two  distributions,  a  large  value  for  D  tends  to  discredit  the  proposed  null 
hypothesis. 

The  computed  values  of  the  test  statistic  D  are  given  in  Tables  C-7  thru  C-9. 
Under  00%  communication  degradation  (see  Table  C-7),  over  all  levels  of  intensity,  the 
computed  K-S  values  show  that  58%  of  the  time  the  null  hypothesis  could  not  be 
rejected.  Table  C-8  shows  that  the  best  results  were  obtained  at  the  15% 
communication  degradation  level  where  9  of  12  cells  passed  the  goodness-of-fit  test  when 
the  method  of  maximum  likelihood  was  used  to  estimate  6.  Finally,  at  the  30%' 
communication  degradation  level,  8  of  the  12  cells  passed  using  maximum  likelihood 
estimates  of  0. 

Over  all  three  levels  of  communication  degradation,  the  null  hypothesis  could  not 
be  rejected  at  the  0.05  level  of  significance  for  24  of  the  36  cells  tested  (67%)  when  F(,(x) 
was  computed  using  maximum  likelihood  estimates  of  0. 

^  Lilliefon,  H.  W.,  ’On  the  Kotmogorov-Sminiov  T*«t  for  the  E^ponentinl  Distnbution  with  Mean  Unknown,’  Journal  of  the 
American  Statirtical  Aaeociation.  Vol  61,  1S69 
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Based  on  the  above  results,  the  specialized  K-S  goodness-of-fit  test  appears  to  have 
qualified  the  exponential  model  as  an  appropriate  fit  to  the  FIST  HQ  service  time  data 
for  fire  request  messages. 

Unlike  the  K-S  test,  the  chi-square  test’s  measure  of  incompatibility  is  based  simply 
on  the  vertical  deviations  between  the  observed  frequencies  and  the  associated  set  of 
expected,  or  theoretical,  frequencies  derived  under  the  assumption  that  the  null 
hypothesis  is  true.  Strictly  speaking,  the  chi-square  test  is  not  appropriate  as  a 
goodness-of-fit  test  for  the  FIST  HQ  artillery  target  intelligence  message  service  time 
data  since  the  exponential  distribution  is  a  continuous  rather  than  a  discrete 
distribution;  however,  it  has  been  justified  by  considering  the  ATI  service  time  data  as 
grouped  into  a  finite  number  of  mutually  exclusive  and  exhaustive  intervals.  The 
observations  in  each  eight-hour  cell  were  grouped  into  r,  arbitrarily  chosen,non- 
overlapping  numerical  categories.  The  data  was  analyzed  in  the  form  of  count  data, 
where  each  count  represents  the  number  of  observations  classified  in  each  of  the  r 
categories. 


The  r  categories  were  chosen  based  upon  a  frequently  used  rule  concerning  the  size 
of  the  expected  frequencies.  The  rule  states  that  no  more  than  20%  of  the  expected 
frequencies  should  be  less  than  5.0,  and  none  can  be  less  than  1.0.  Based  on  this  rule, 
the  low,  medium  and  high  intensity  levels  for  ATI  messages  were  grouped  into  8,  7,  and 
4  categories,  respectively. 


The  null  hypothesis  to  be  tested  is  that  the  probability  density  function  of  the 
observed  random  variable,  ATI  message  service  time,  is  some  fo(x),  set  against  the 
alternative  hypothesis  that  the  probability  density  function  of  the  observed  random 
variable  is  different  from  f0(x). 


The  test  statistic  computed  to  make  the  aforementioned  comparison  is  denoted  as 
X^.  This  statistic  may  be  defined  as 


(8) 


where  fo  =  an  observed  frequency,  and  fj  =  an  expected,  or  theoretical,  frequency  for  a 
particular  category.  Associated  with  this  statistic  is  a  number  of  degrees  of  freedom 
which  is  equal  to  the  number  of  categories  for  each  communication  degradation  level 
minus  one.  In  addition,  for  this  analysb,  the  single  parameter  of  the  exponential 
distribution  that  was  estimated  from  the  data  prior  to  determining  the  expected 
frequencies  also  had  to  be  subtracted  off  the  total  number  of  degrees  of  freedom. 


From  Eq.  (8)  it  can  be  seen  that  if  every  observed  frequency  is  exactly  equal  to  the 
corresponding  expected  frequency,  then  the  value  of  x^  is  0.0  (the  smallest  possible  value 
of  x*).  The  larger  the  discrepancies  between  the  observed  and  expected  frequencies,  the 
larger  the  value  of  x^-  The  null  hypothesis  was  tested  at  the  0.05  level  of  significance. 
This  means  that  if  the  null  hypothesis  were  true,  the  probability  of  observing  a 
computed  value  greater  than  the  critical  (tabulated)  x*  value  would  be  0.05.  The 
computed  x*  values  are  given  in  the  body  of  Table  C-12.  The  appropriate  critical  x* 
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values  are  12.592,  11.070,  and  5.991  for  6,  5,  and  2  degrees  of  freedom,  respectively. 
Approximately  78%  (7  of  9  cells)  of  the  time  the  null  hypothesis  could  not  be  rejected 
when  maximum  likelihood  estimates  of  0  were  used  in  computing  the  expected 
frequencies  needed  for  evaluating  the  test  statistic.  The  two  conditions  under  which  the 
X*  statistic  was  significant  were  "perfect”  communications  (00%  communication 
degradation)  with  low  intensity  and  15%  communication  degradation  with  medium 
intensity.  These  large  computed  x^  values  do  not  appear  to  be  due  to  anomalies  within 
the  observed  data,  but  rather  to  the  unfortunate  instability  of  the  chi-square  statistic 
since  its  value  is  affected  by  the  number  of  categories  within  a  cell  and  the  width  of 
each  category. 

The  hypothesis  that  the  population  probability  distribution  of  ATI  message  service 
times  is  exponential  was  tested  and,  overall,  could  not  be  rejected.  Thus,  at  the  0.05 
level  of  significance,  the  exponential  distribution  can  reasonably  be  used  as  a  model  for 
FIST  HQ  service  time  for  artillery  target  intelligence  messages. 

m.  MODELING  RESULTS 

Modeling  efforts  of  the  FIST  HQ  fire  request  and  artillery  target  intelligence 
message  service  time  data  collected  during  the  first  CPX  Research  Facility  experiment 
indicate  that  the  distribution  of  the  service  times  of  these  message  types  can  be  initially 
fit  with  a  three-parameter  Weibull  distribution.  Under  the  assumptions  that  one,  the 
shape  parameter,  of  the  distribution  is  equal  (or  close)  to  1.000,  and  two,  the  location 
parameter,  u,  takes  on  the  value  of  the  smallest  order  statistic  of  each  of  the  data  sets 
modeled,  the  distribution  of  FIST  HQ  service  times  for  FR  and  ATI  messages  is  easily 
transformed  into  a  single  parameter  exponential  distribution.  The  maximum  likelihood 
estimate  o^  the  parameter  0  of  thb  distribution  is  easily  computed  as 

0  =  ( 1.0/  ( Xi/N)  j,  where  the  y-,  are  the  adjusted  service  times  and  N  is  the  sample 
i-l 

size  of  the  data  set. 

There  are  certain  trade-offs  made  by  electing  to  use  the  exponential  distribution 
rather  than  the  three-parameter  Weibull  distribution  to  model  FR  and  ATI  message 
service  times.  The  Weibull  distribution  performed  well  in  describing  the  rise-and-fall 
pattern  exhibited  by  many  of  the  data  sets,  accounting  for  both  increasing  and 
decreasing  service  times.  However,  it  is  considered  a  "short  tailed”  distribution  and,  as 
such,  did  not  do  well  in  modeling  extreme  value  observations.  On  the  other  hand,  the 
exponential  distribution  is  strictly  a  monotonically  decreasing  function  and  assumes  a 
constant  service  time  rate.  Although  the  exponential  density  function  fell  off  rapidly  as 
service  time  increased,  in  most  cases,  the  steepness  of  the  fall  was  very  gradual  as 
evidenced  by  the  small  values  of  0  (all  <  1.000)  that  were  computed.  This  long  tail  was 
able  to  accommodate  the  large  service  times  observed  in  some  cells. 

A  specialized  Kolmogorov-Smirnov  test  and  the  chi-square  good ness-of- fit  test  were 
used  to  compare  the  distribution  of  the  transformed  service  times  with  the  theoretical, 
or  hypothesized,  distribution  (i.e.,  the  single  parameter  exponential  distribution)  for  fire 


request  and  ATI  messages,  respectively.  For  each  of  the  cells  modeled,  the  unknown 
parameter  0  had  to  be  estimated  from  the  data  before  the  tests  could  be  performed. 
Results  showed  that  when  utilizing  maximum  likelihood  estimates  of  B,  67%  of  the  fire 
request  message  cells  (two-hour)  could  not  be  rejected  at  the  0.05  level  of  significance, 
while  7S%  of  the  ATI  message  cells  (four-hour)  could  not  be  rejected  at  the  same  level 
of  significance. 

While  a  classical  parametric  model  provided  a  reasonable  fit  to  the  distribution  of 
the  observed  FIST  HQ  service  times  for  FRs  and  ATIs,  it  should  also  be  realized  that  it 
does  obscure  situations  that  can  arise  at  the  FIST  HQ  and  ultimately  influence  its 
message  service  time.  Problems,  such  as:  the  prioritization  of  particular  types  of 
missions  (like  Copperhead),  the  limited  queue  size  of  the  FIST  DMD  (it  can 
accommodate  only  16  messages  at  any  one  time),  the  fact  that  messages  may  be  turned 
away,  and  a  varying  message  arrival  rate  at  the  FIST  HQ  (due  to  a  change  in  the  level 
of  intensity  of  the  battle)  could  all  be  handled  by  the  development  of  an  efficient 
queueing  system  model  for  the  FIST  HQ  service  time. 

rv.  CONCLUSIONS 

A  discussion  of  the  procedure  by  which  the  FIST  HQ  service  time  data  for  FRs  and 
ATIs  was  modeled  and  the  techniques  utilized  to  test  the  fit  of  the  hypothesized  model 
has  been  provided.  It  was  found  that  while  the  FIST  HQ  was  operating  under  the 
forward  observer  (FO)  review  control  mode,  the  service  times  for  fire  request  messages 
and  artillery  target  intelligence  messages  were  reasonably  modeled  by  a  single  parameter 
exponential  distribution. 

The  exponential  distribution  generated  from  the  data  is  easy  to  work  with  and  it 
has  shown  its  capability  of  accounting  for  differences  among  levels  of  intensity  and 
communication  degradation,  team,  and  replicate.  However,  it  does  "mask”  other 
conditions  existent  at  the  FIST  HQ  which  could  also  significantly  influence  message 
service  time.  Despite  this  fact,  the  assumption  that  message  service  times  at  the  FIST 
HQ  are  exponentially  distributed  can  be  successfully  used  as  input  into  the  development 
of  a  queueing  system  model  that  would  more  thoroughly  represent  the  conditions  that 
influence  the  FIST  HQ  message  service  time. 
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TABLE  C-IA. 


POINT  E^MATES  POR  THE 
3-PARAMEIER  ^BULL  DISTRTBUnON 
FOR  FIRE  REQUEST  MESSAC^ 
(00%  Degradatbn) 


Panuneters 

R^ 

Team 

Intensity 

a 

V 

L 

4.027 

1.536 

10.667 

3 

M 

4.708 

1.176 

10.889 

1 

H 

12.925 

1.387 

9.200 

L 

4.274 

1.620 

8.500 

4 

M 

6.969 

1.361 

4.000 

H 

8.197 

1.203 

4.000 

L 

10.345 

2.040 

7.500 

3 

M 

6.335 

2.079 

2 

H 

14.444 

2.128 

L 

4.396 

2.157 

PSl 

4 

M 

5.319 

1.016 

8.000 

H 

9.685 

1.919 

TABLE  C-IB. 


MAXIMUM  LIKELIHOOD  ESTIMATES  FOR  THE 
3-PARAMEIER  WEIBULL  DISTRIBUTION 
FOR  FIRE  REQUEST  MESSAGES 
(0(M  Degradation) 


Rq> 

Team 

Intensity 

Parameters 

a 

& 

V 

L 

8.610 

6.447 

11.000 

3 

M 

11.263 

3.190 

11.000 

1 

H 

8.488 

2.450 

11.000 

1 

L 

5.414 

1.610 

9.000 

4 

M 

1.050 

1.175 

5.000 

H 

1.992 

1.000 

5.000 

L 

10.020 

1.638 

9.000 

3 

M 

7.299 

1.926 

10.000 

H 

7.281 

1.830 

8.000 

A 

L 

3.399 

1.691 

4 

M 

3.599 

1.081 

8.000 

H 

6.903 

1.717 
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TABLE  C-2A. 


POIKT  ESTIMATES  FOR  THE 
3-PARAMEIER  WHBULL  DISTRIBUTION 
FOR  FIRE  REQUEST  MESSAGES 
(1S%  Degradation) 


TABLE  C-2B. 


TABLE  C-3A. 


POINT  ESmMATES  FOR  THE 
3-PARAMEIER  WEIBULL  DISTRIBUnON 
FOR  FIRE  REQUEST  MESSAGES 
(30%  Degredatbn) 


Rep 

Team 

Intensity 

Parameters 

a 

V 

L 

41.781 

4.346 

25.000 

3 

M 

9J84 

0.787 

11.923 

1 

H 

13.812 

.995 

10.500 

L 

11.801 

1.141 

4.333 

4 

M 

13.043 

2.187 

1.000 

H 

14.205 

.976 

3.500 

L 

21.127 

.871 

16.000 

3 

M 

23.018 

1.809 

5.000 

2 

H 

22.601 

1.299 

6.875 

L 

7.627 

2.049 

3.667 

4 

M 

5.000 

1.410 

6.500 

H 

8.632 

.850 

11.826 
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TABLE  C-4. 


MAXIMUM  LTKELTHCX)DESrnMATES  OF  0  FOR  THE 
SINGLE  PARAMETER  EXPONENTIAL  DISTRIBUTION 
FOR  FIRE  REQUEST  MESSAOTS 
(00%  Degradation) 


TABLE  C-6. 


maximum  LIKELIHOOD  ESTIMATES  OF  »  FOR  THE 
SINGLE  PARAMETER  EXPONENDAL  DISTRIBUTION 
FOR  FIRE  REQUEST  MESSAGES 
(30%  Degradation) 


Rq) 

Team 

Parameter 

Intcn&ity 

9 

L 

.046 

3 

M 

H 

.072 

1 

L 

.115 

4 

M 

.128 

H 

.107 

L 

.147 

3 

M 

.060 

H 

.059 

2 

■m 

M 

H 


TABLE  C-7. 


COMPUTED  KOLMOGOROV-SMIRNOV  TEST  STATISTIC  VALUES 
USING  MAXIMUM  LIKELIHOOD  ESTIMATES  FOR  THE 
SINGLE  PARAMETER  EXPONENTIAL  DISTRIBUTION 
FOR  HRE  REQUEST  MESSAGES 
(00%  Degradation) 


TABLE  C-8. 


COMPUTED  KOLMOGOROV-SMIRNOV  TEST  CTATISTIC  VALUES 
USING  MAXIMUM  LIKELIHOOD  ESTIMATES  FOR  THE 
SINGLE  PARAMETER  EXPONENTIAL  DISTRIBUTION 
FOR  iTRE  REQUEST  MESSAGES 
(15%  Degradation) 


R^ 

Team 

Intensity 

Computed  K>S  Values  (D) 

L 

J43 

3 

M 

20S 

1 

H 

2U 

L 

.400 

4 

M 

.369* 

H 

.235 

L 

.411 

3 

M 

.292* 

2 

H 

.202 

L 

.324 

4 

M 

.416* 

H 

.113 

•Indicates  the  value  of  D  is  significant  at  a  level  of 
significance  of  .05. 
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TABLE  C-9. 


COMPUTED  KOLMOGOROV-SMIRNOV  TEST  STATICTIC  VALUES 
USING  MAXIMUM  LIKELIHOOD  ESTHMATES  FOR  THE 
SINGLE  PARAMETER  EXPONENTIAL  DISTTRIBUTION 
FOR  HRE  REQUEST  MESSAGES 
(30%  D^radatbn) 


^Indicates  the  value  of  D  is  significant  at  a  level  of 
significance  of  .05. 
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TABLE  C-IOA. 

POTM*  ESTIMATES  FOR  THE  S-PARAMETER 
WHBULL  DISTRIBUTION  FOR  ATIs 
(D^radation  by  IntenslQr) 


TABLE  C-11. 


MAXIMUM  LIKELIHOOD  ESTIMATES  FOR  THE 
SINGLE  P  ARAMEIER  EXPONEND  AL  DISTRIBUTION 

FOR  ATIs 

(D^radation  by  IntensiQr) 


Communication 
Degradation  (%) 

IntensiQr 

Parameter 

0 

L 

.150 

00 

M 

.194 

H 

.240 

L 

.180 

15 

M 

.096 

H 

.112 

L 

.104 

30 

M 

.103 

H 

.041 

vv 


TABLE  012. 


COMPUTED  ^  VALUES  USING 
MAXIMUM  LIKELIHOOD  ESTIMATES  FOR  THE 
SINGLE  PARAMETER  EXPONENDAL  DISTRIBUTION 

FOR  ATIs 


Communication 
Degradation  (%) 

Intensi^ 

Computed  Values 

Degrees  of 
Freedom 

L 

24.354* 

6 

00 

M 

7.672 

S 

H 

3.882 

2 

L 

6 

15 

M 

5 

H 

.977 

2 

L 

11.026 

6 

30 

M 

10.731 

5 

H 

2.998 

2 

•Indicates  a  significant  value  at  the  .05  level  of  significance. 
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APPENDIX  D. 


RELATIONSHIP  BETWEEN  THE  WEIBULL  AND  EXPONENTIAL  DISTRIBUTIONS 


Suppose  that  the  random  variable  X,  which  represents  the  FIST  HQ  service  time, 
has  a  three-parameter  Weibull  distribution,  then  the  probability  density  function, 
fx  (x)  of  X,  has  the  following  form; 


a  a  a 


if  X  >  i/ 


fx  (x)  = 


if  X  <  I/;  0,ot>  0;  (Dl) 

1/  >  0. 


The  three  constants  >  0,  a  >  0,  and  */  >  0  are  the  parameters  of  the 
distribution.  The  parameter  P  determines  the  shape  of  the  density  function.  The 
parameter  a  is  a .  scale  parameter  specifying  the  100  (1  -  exp  (-l)]th  distribution 
percentile  of  X  -  i/.  The  parameter  i/  will  assume  the  smallest  possible  value  of  the 
random  variable  X.  It  may  be  thought  of  as  a  location  (or  threshold)  parameter  where  a 
message  will  be  serviced  before  time  i/  with  probability  0. 

Weibull  random  variables  may  be  easily  transformed  to  exponentially  distributed 
random  variables.  An  exponential  model  can  be  utilized  if  an  appropriate 
transformation  of  the  FIST  HQ  service  time  data  is  made.  Suppose  ordered  service 
times  X(i),  •  •  •  ,  Xj^j  from  a  sample  of  size  n  are  observed  and  the  observations  are  from 
a  population  of  Weibull  random  variables  with  known  shape  and  location  parameters 
0  and  i^,  respectively,  and  unknown  scale  parameter  a.  Then  y^ij  =  (X||)  -  i/)^,  ■  •  , 
y|^j  =:  will  be  considered  ordered  o^ervations  from  an  exponential 

distribution  with  unknown  scale  parameter  ^  =  (-^)  . 

Q 

This  transformation  is  based  on  the  assumption  that  0  and  v  are  assumed  to  be 
known  for  a  given  set  of  treatment  combinations.  Based  on  the  maximum  likelihood 
estimates  for  0  and  u  (see  Tables  C-lB,  C-2B,  C-3B  and  C-lOB),  the  value  for  0  will  be 
set  equal  to  1.000,  while  u  will  be  set  equal  to  the  smallest  order  statistic  of  each  data 
set.  Thus,  the  unknown  scale  parameter  of  the  exponential  distribution  becomes 


To  verify  the  assertion  that  the  Weibull  distribution  is  equivalent  to  the 

exponential  distribution  with  parameter  9  =  —  under  the  above  assumptions,  the 

.®  .  .  . 

ordered  observations  from  the  exponential  distribution  yjjj  =  (xjij  - 1/)",  ..., 
y(r)  =  (X{r)  -  can  *^6  rewritten  as  y^,)  =  (x,,)  -  i/),  ...,  =  (X(,)  -  v). 

Thus  Y  =  ( (X  -u)  I  =g(X)  is  a  general  transformation  of  the  continuous  random 
variable  X;  we  would  like  to  obtain  the  distribution  of  Y  from  the  distribution  of  X. 
Since  the  observations  y^jj,  ...,  y^^j  are  from  the  exponential  distribution,  then  Y  =  g(X) 
is  a  strictly  decreasing,  differentiable  transformation.  A  strictly  decreasing 


transformation  of  the  form  y  =  g(x)  implies  that  if  g(x,)  >  glxo)  then  <  Xj. 
Differentiable  implies  (d/dx)g(x)  exists  for  every  x.  If  y  =  g(x)  =  (x  -  i/)  is  strictly 
decreasing  and  differentiable,  there  exists  a  decreasing,  differentiable  function  h(y)  such 
that  X  =  h(y)  =  y  +  i/. 

Thus,  the  density  function  of  Y  =  g(X)  becomes, 

Fy  (y)  =  p  {Y  <  y}  =  P  {g  (X)  <  y}  =  P  {X  >  h(y)  }  =  1  -  Fx  (h(y)) .  (D2) 


Substituting  g(X)  =  (X  -  i/)  and  h(y)  =  y  +  i/,  equation  D2  becomes, 

Fy  (y)  =  P  {Y  <  y}  =  P  {(X-  y)  <  y}  =  P  {X  >  (y  +  y) }  (D3) 

=  1  -  Fx  (y  +  t/) . 

If  both  sides  of  the  equality  are  differentiated  with  respect  to  y,  the  probability  density 
function  of  Y  is 

fy  (y)  =  ( -  fx  (My)))  =  1  1  tx  (h(y)) ,  (D4) 

since  (d/dy)  h(y)  <  0. 


(D5) 


Substituting  for  h(y)  gives 

f,  ,y I  =  (  -  f,  ,y  +  =  1  1  ,,,  (y  +  .)  . 

Hence,  the  probability  density  function  of  Y  is  ^ 

|],D5, 

Recall  that  ^  =  1.000,  then 


thus 


frly)  =  :^«p|-(i)  ). 

It  was  previously  stated  that  =  — ,  therefore,  substituting  in  equation  D8  gives 


(D8) 


fy  (y)  =  0  exp  (  -  (^)  ) . 


(D9) 


Therefore,  the  distributional  form  of  the  transformed  random  variables  y|j),  ...,  yj,j  is 
exponential  with  unknown  scale  parameter  9,  while  the  original  observations  were 
considered  to  be  from  a  population  of  Weibull  random  variables. 
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